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i> ■ ABSTRACT 

High-resolution N-body simulations of hierarchical clustering in a wide variety of 
O . cosmogonies show that the density profiles of dark matter halos are universal, with 

low mass halos being denser than their more massive counterparts. This mass-density 



in 



correlation is interpreted as reflecting the earlier typical formation time of less massive 
objects. We investigate this hypothesis in the light of formation times defined as the 
epoch at which halos experience their last major merger. Such halo formation times 
f"^ ' are calculated by means of a modification of the extended Press &: Schechter formalism 

OO ■ which includes a phenomenological frontier. Am, between tiny and notable relative mass 

^nq ■ captures leading to the distinction between merger and accretion. For Am ~ 0.6, we 

confirm that the characteristic density of halos is essentially proportional to the mean 
Q>^ , density of the universe at their time of formation. Yet, proportionality with respect to 

the critical density yields slightly better results for open universes. In addition, we find 
Qh' that the scale radius of halos is also essentially proportional to their virial radius at the 

O ' time of formation. 

5— ( ' 

-*— > ' We show that these two relations are consistent with the following simple scenario. 

Ci • Violent relaxation caused by mergers rearranges the structure of halos leading to the 

same density profile with universal values of the dimensionless characteristic density 
k> ] and scale radius. Between mergers, halos grow gradually through the accretion of 

H ' surrounding layers by keeping their central parts steady and expanding their virial 

radius as the critical density of the universe diminishes. 

Subject headings: cosmology: theory - galaxies: halos ~ galaxies: formation 



1. Introduction 

High-resolution N-body simulations of hierarchical clustering in the standard CDM cosmogony 
carried out by Navarro, Frenk, & White (1996) show that the spherically averaged equilibrium 
density profiles of dark matter halos with masses ranging from dwarf galaxy to rich cluster scales 
are well fitted (see however Moore et al. 1997) by the formula 

p{x) x^ 

Pcvit X{x + Xs)'' 



In equation (|l]), pcrit is the critical density of the universe, x = r/R is the radius scaled to the 
so-called virial radius R, and 5c = Pel Pent and Xs = r^jR are two dimensionless parameters giving 
the characteristic density and scale radius, respectively, of the density profile. These latter two 
parameters are linked through the relation 

. ^ 200 x^ 

' 3 [ln(l+Xs-^)-(l + x.)-i]' ^ ^ 

arising from the steadiness condition that the mean density within R is equal to 200 x Pcrit- Thus, 
the dimensional density profile of a halo with mass M at a given time t (the latter two quantities 
fixing the values of R and pcrit in a given cosmogony) is governed by one single free parameter. 
Note that the inverse of Xg is a direct measure of the halo concentration. 

More recently, Navarro, Frenk, & White (1997, hereafter NFW), and in independent work, Cole 
&: Lacey (1997) and Tormen, Bouchet, & White (1997), have shown that expression (|I|) provides 
equally good fits to the density profile of dark halos in a number of other cosmogonies, including 
flat and open models, with or without a cosmological constant, and with different initial power 
spectra of Gaussian density fluctuations. In all the cosmogonies investigated the parameter b^ has 
been found to correlate with mass in such a way that low mass halos are denser than those of high 
mass. This mass-density correlation is interpreted as reflecting the earlier typical formation time 
of less massive objects. As shown by NFW, the correlation is well described by a simple model in 
which the characteristic density pc of a halo of present mass Mq is proportional to the mean density 
of the universe at the corresponding formation redshift Zf , or equivalently, 

5e = CJ]o[l + 2f(Mo)]'. (3) 

To compute Zf (Mq), NFW used the expression derived by Lacey & Cole (1993, LC) in the framework 
of the Press & Schechter (1974, PS) prescription for the cumulative probability that the mass of a 
halo following single M(f) tracks reduces to some fraction of its present mass, /, taken by NFW 
as a free parameter. They find that for / < 0.01 the predicted typical mass-density relations fit all 
their simulations reasonably well for essentially the same proportionality factor C . Although this 
result strictly refers to present-day halos, it should also apply to halos at any redshift for scale- free 
cosmogonies and those in which the evolution of structure is close to being self-similar. 

In spite of this remarkable result one cannot overlook the fact that the distribution of formation 
times based on single M(t) tracks is not fully adequate for estimating the time at which a parent 
halo reaches, for the first time, a fraction / of its present mass (cf. LC). Moreover, the fact that / 
must be less than or equal to 0.01 poses two problems. Firstly, it leads to an ambiguous definition 
of the formation time, since a progenitor with M < O.SMq is not necessarily along the main lineage. 
Second, it is difficult to understand how the present structure of a halo can bear any relationship to 
the epoch in which some progenitor reached such a small fraction of the current halo mass. More 
importantly, the definition of the formation time in the LC clustering model does not distinguish 
between notable mass increases occurring more or less abruptly in time. Major deviations from 



equilibrium and subsequent violent relaxation take place only when halos of comparable masses 
merge, while tiny mass captures have a negligible effect on the structure of the capturing systems. 
Numerical simulations of hierarchical clustering (Cole & Lacey 1997; Thomas et al. 1997) indeed 
show that halos with no evidence of a recent major merger are in steady state within R, despite 
the fact that they are continually capturing small halos. 

Kitayama & Suto (1996) have attempted to describe the formation and destruction of halos 
within the extended PS prescription by differentiating between notable and tiny relative mass 
captures. Their model lacks, however, a consistent definition for the formation of halos because all 
halo captures involved in the same common merger are counted separately as giving rise to different 
new halos. A similar, but fully consistent approach, has been followed independently by Manrique 
& Salvador-Sole (1995, 1996, hereafter MS95 and MS96). These authors have developed a semi- 
analytical clustering model within the framework of the peak formalism, hereafter referred to as the 
CUSP (Confluent System of Peak trajectories) model, which distinguishes naturally between major 
and minor mergers, hereafter simply referred to as (true) mergers and accretion. This allows one 
to define unambiguously the halo formation and destruction times corresponding, respectively, to 
their last and next merger. Unfortunately, to be fully satisfactory the CUSP model requires a more 
accurate expression for the peak-peak correlation at small separations than is presently available 
(Manrique et al. 1997, M&CO). 

In the present paper, we develop a self-consistent modification of the LC model which, drawing 
inspiration from the CUSP model, differentiates merger from accretion. This model, which retains 
the simplicity and good predictive properties of the original model (Lacey &: Cole 1994) while in- 
cluding better motivated formation and destruction time estimates, is used to investigate the origin 
of the empirical mass-density and related mass-radius correlations, as well as their implications for 
the evolution of halo structure in hierarchical cosmogonies. The paper is organized as follows. The 
modified LC model is presented in § ^ It is applied in § ^ to the study of the empirical mass-density 
and mass-radius correlations. The results of this analysis are summarized in § 0. 



2. Merger vs. Accretion and the PS Formalism 

A halo survives as long as it evolves by accretion, or equivalently, as long as it captures only 
relatively tiny systems. Otherwise, it merges, which automatically leads to its destruction. Note 
that when a halo is captured by one that is more massive it merges and is destroyed in the event, 
but the capturing halo may survive provided that the captured mass is relatively small. Only those 
events in which all the initial halos merge and are destroyed give rise to the formation of new halos. 

The preceding definitions do not affect the abundance of halos at a given time, only the 
description of their growth. It is, therefore, not surprising that the CUSP model also distinguishing 
between merger and accretion predicts a halo mass function (MS95) that is highly similar to the 
PS one as in the LC model. Accordingly, in the modified LC clustering model we propose, the mass 



function is equal to the PS mass function 



NiM,t) = Ni^ciM,t) 
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(4) 



where po is the present mean mass density of the universe, (5coii(t) is the critical overdensity for 
collapse at t linearly extrapolated to the present time, and a{M^ is the current r.m.s. overdensity 
on spheres encompassing a mass M. 

The instantaneous merger rate for halos of mass M at t per infinitesimal range of final masses 
M' > M, or specific merger rate, predicted by the CUSP model (MS96) is also close to the 
corresponding rate predicted by the original LC model, down to some value Am of the relative 
captured mass AM/M = (M' — M)/M, where it shows a sharp cutoff. This cutoff reflects the fact 
that, in the CUSP model, captures of small mass halos relative to M are not computed as mergers, 
but simply contribute to accretion. In contrast, the specific merger rate predicted by the LC model 
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keeps on increasing monotonically for small AM/M because any mass capture is regarded, in this 
model, as a merger, and the number density of small mass halos diverges. Following this result we 
modify the original LC model by including a threshold Am in the relative mass captured by a halo 
for such an event to be considered a merger, smaller mass captures only contributing to continuous 
accretion. With this modification the new specific merger rate takes the form 



r'^iM ^ M',t) 
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if M < M' < M(Am + 1) 
M',t) if M(Am + 1) < M', 



(6) 



while the total mass increase rate for halos of mass M at t, rmass(-^;i) = dM/dt, splits into two 
contributions, one arising from mergers, or mass merger rate, 



,{M,t) 



and the other arising from accretion, or mass accretion rate, 



AMr™(M^M',t)dM', 
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AM rfc(M -^ M',t) dM'. 



(7) 



(8) 



As shown by M&CO, the mass function, the specific merger rate, and the mass accretion rate 
determine the behavior of the entire CUSP model. This is also the case for the modified LC model. 



The specific merger rate determines the mass merger rate (eq. [^), as well as the destruction rate 
of halos with mass M at t 

POO 

r'^{M, t) = / r'^iM -^ M', t) dW . (9) 

Likewise, the formation rate can be written as 

r^[M{t),t] = dlniV[M(t),t] + r^^lMit^] + dMr^^^.{M,t)\^^^^,^, (10) 

from the conservation equation for the number density of halos per unit mass along mean mass 
accretion tracks, M{t), solution of the differential equation 

^=<ass[M(t),t]. (11) 

Finally, the distributions of formation and destruction times in the modified LC model are given 
by expressions identical to those in the CUSP model (see M&CO). In particular, the distribution 
of formation times for halos at to with masses between Mq and Mq + 6Mo, with 6Mq arbitrarily 
small, takes the form 

*'<*> ^1^)^ = --'[^'W-*! "^{-C '•'[''('''•''1 '''}• <i^> 

with M(t) the mass of these halos at t calculated along their mean mass accretion tracks. The 
median of this distribution is adopted as the typical halo formation time. 

Before concluding this section, we should clarify the fact that the distinction adopted between 
merger and accretion is not motivated by the results of A^-body simulations, but obeys the desire to 
differentiate schematically the dynamic effects on halo structure of tiny and notable relative mass 
captures. Note also that while the merger cutoff in the CUSP model arises naturally from the 
peak ansatz and the assumed distinction between merger and accretion (see MS96), the threshold 
for merger. Am, in the present modified LC model is a free phenomenological parameter which, for 
simplicity, will be considered independent of M and t (one assumption implies the other in scale- 
free universes). Strictly speaking, the rearrangement of a halo after the merger of two progenitors 
depends on the relative gain of energy per unit mass rather than simply on the relative mass 
increase. However, as the former quantity is largely determined by the latter, this simplifying 
assumption is justified. 



3. The (5c (Mo) and x^iMo) Correlations 

Next we investigate the possible origin of the mass-density and mass-radius correlations ex- 
hibited by halos in hierarchical universes in the light of the modified LC model developed in the 



preceding section. To do this we use the numerical data of NFW, which comprises the eight dif- 
ferent cosmogonies hsted in Table 1. The first column of this Table lists the power spectra, while 
columns (2) and (3) list the values of Qq and Aq = A/(3Hq), respectively. We list in column (4) the 
present r.m.s. density fluctuation in 8h~^ Mpc spheres, ag. To facilitate the comparison among the 
cosmogonies, masses are scaled to the values of the present characteristic mass M^, defined through 
a{M^:) = (^coii(*o)) which are listed in column (5) of Table 1. All the models have h = 0.5, except the 
ACDM model which has h = 0.75, with the Hubble constant defined as Hq = 100/ikms~^Mpc~^. 

In Figure 1 we show the best fits (by a standard least squares minimization in logarithmic 
units) to the empirical 6c{Mq) correlation obtained using the fitting formula (y) for three different 
values of Am. The value 0.6 corresponds to the best overall fit when Am is varied from 0.1 to 0.9 
in steps of one tenth. This value is also favored individually by the three flat power-law spectrum 
models with n = —1, —0.5, and 0, which are those that best discriminate among the different 
predictions. The two open scale-free models favor a value of Am = 0.5, while the remaining power 
spectrum model and the two CDM models favor Am = 0.7, although marginally. In other words, 
as it is apparent from Figure 1, a value of Am ~ 0.6 gives reasonably good fits in each of the 
cosmogonies investigated. In contrast, the predictions corresponding to the extreme values 0.1 and 
0.9 do not describe the numerical data well in practically any case. 



Given the formation time distribution function (eq. |12]) and relation @ we can readily 
derive the distribution functions of log(5c) for any value of Mq. Figure 2 shows the distributions 
obtained in the SCDM model for five different values of Mq, other cosmogonies giving qualitatively 
similar results. They are in good overall agreement with the empirical distributions of points: the 
maxima are near to the values of log((5c) corresponding to the median formation redshifts, and 
the spreads have the right magnitude and show a trend to diminish with increasing Mq. This 
indicates that relation (|3|) also applies to individual halos and that their characteristic density, pc, 
remains essentially equal to C times the mean cosmic density at their time of formation. Note that, 
according to the PS mass function, low mass halos are severely underrepresented in the empirical 
samples with respect to more massive ones, indicating that the selection of the former has been 
much stricter. In this manner, earlier formation times may have been artificially favored, since 
the older the halos, the better they satisfy the requirement of having a relaxed appearance. This 
might explain the slightly smaller scatter shown by the empirical distributions for small mass halos. 
This effect and the slight bias also introduced by our simple fitting procedure (we have assumed 
constant, symmetrical errors) might affect to some extent the quantitative results of the fits, but 
the previous conclusions should prevail. 

The values of C listed in column (6) of Table 1 show a much wider variation with the cosmogony 
(an overall factor 100) than in NFW (only a factor two; see their Table 1). Although the possible 
biases mentioned above might in part be responsible for this variation, the marked departure from 
a hypothetical common value shown by the values of C in open cosmogonies seems real. We 
have investigated the possibility of reducing the scatter in the Oq < 1 cases by devising a slightly 
different model which has the added value of providing a straightforward physical interpretation of 



the empirical mass-density correlation. In the new model, the characteristic density, pc, of halos 
with current mass Mq is assumed to be proportional to the critical density of the universe at their 
time of formation, instead of to the mean cosmic density at that time. With such a proportionality, 
not only does pc remain fixed from the time of halo formation, but also the initial value of 5c is 
universal (i.e., independent of mass and time in self-similar universes). From the form this fitting 
model adopts in dimensionless units 

it is apparent that the value of 6c when halos form is equal to the proportionality factor 6ci- 

The best overall fit of the empirical data with the model (|l^) is obtained again for Am = 0.6. 
As can be seen from Figure 1, the fits in the open cases are slightly better than in the original 
model (^), while the two models give, of course, identical results in the f^o = 1 cases. The overall 
variation shown by the proportionality factor 6c{ in different cosmogonies has diminished, although 
a trend of 6ci with cosmogony is still present. We note that some theoretical studies predict a 
dependence on the cosmogony of the typical halo density profiles resulting from violent relaxation 
(e.g., Syer k White 1997). 

Relation (|2|) between the dimensionless parameters 6c and Xs tells us that the value of x^ 
shown by halos at their time of formation, hereafter denoted by Xgf, is also universal. The values 
of Xsf inferred from those of 5cf drawn from the previous fits are listed in column (8) of Table 
1. The universality of Xgf is equivalent to stating that the dimensional scale radius rg of halos 
at their time of formation is proportional, with universal proportionality factor equal to Xgf, to 
their virial radius R at that epoch. This raises the question: is the scale radius rg of current halos 
also proportional, with identical proportionality factor, to their virial radius R at their time of 
formation? Or equivalently, does the value of rg for current halos of mass Mq coincide with the 
value this parameter had when they formed, as for pc? The answer to these questions is not trivial 
since it depends on the relation between the initial and current halo masses, that is, on the typical 
mass accreted since their formation. 

The modified LC model allows us to estimate the mass accreted by halos. Hence, we can 
readily check the previous proportionality, which in dimensionless units takes the form 

R{M[zdMo),Mo]} 
^^"^^' i?(Mo) ■ ^ ^ 

In Figure 3 we show the results of directly fitting this model to the Xs{Mo) empirical correlation. 
Once again we find the best overall fit for Am = 0.6. More importantly, the best values of Xgf, 
listed in column (9) of Table 1, are in fairly good agreement (just slightly larger on average) with 
those listed in column (8). Note that, despite the convoluted calculations involved (given a halo of 
current mass Mq one must first calculate its formation redshift, then, using the modified LC model, 
its mass at Zf, and finally, through the steadiness condition, the corresponding value of R), the fits 
are as good as those obtained for the mass-density correlation. 



The agreement, case by case, between the two independent values of Xgf given by the corre- 
lations Sc{Mq) (through the relation Q) and Xs{Mq) (through our clustering model) supports the 
overall validity of the modified LC model for A^ ~ 0.6, and of the theoretical relations (13) and 



(14). The physical interpretation of the latter can be summarized as follows: i) the values of the 
dimensionless parameters, 6ci and Xgf , characterizing the radially averaged density profiles of halos 
at formation are universal, and ii) the values of the corresponding dimensional parameters, pc and 
Tg, remain fixed as long as halos evolve by accretion. 

The fact that for a given halo the values of parameters pc and rg are set when it forms tells 
us that the only effect of accretion is the gradual expansion of the halo virial radius R in order to 
permanently satisfy the steadiness condition. We have directly tested this corollary by comparing 
the mass increase experienced by halos since their formation predicted by the modified LC model, 
with the mass increase that results from taking halos with a density profile of the form (|^), with 
fixed values of pc and rg, and progressively increasing the virial radius R as pcrit diminishes. As 
expected, we have found a fair degree of agreement between the two mass evolutions, the maximum 
departure being less than 35% in any one cosmogony. 



4. Conclusions 

The 5c{Mq) correlations predicted by the modified LC model for Am = 0.6, assuming the 
relation (^, match the empirical data as well as the NFW predictions for / = 0.01. We therefore 
confirm, with a more compelling formation time estimate, the claim by these authors that the 
characteristic density shown by dark halos in equilibrium is proportional to the mean density of 
the universe at the time they form. We want to stress that while the two different formation time 
estimates give similarly good fits, this does not imply that the difference in their definitions is merely 
formal: for any given halo mass, the typical formation redshifts used by NFW are appreciably larger, 
typically by a factor of two, than those obtained in the model presented here. We have also shown 
that the fits for the open models can be improved further if one assumes instead a proportionality 
with respect to the critical density of the universe at halo formation. 

The modified LC model presented in § ||, together with the definitions of halo formation and 
destruction times with which it deals, relies on the schematic differentiation of the dynamic effects 
of tiny and notable mass captures. According to this scenario, the structure of halos would be fixed 
through violent relaxation in the last major merger which they had experienced, while between two 
consecutive major mergers halos would remain essentially unaltered, mass accretion only producing 
a progressive expansion of their envelope as new surrounding layers fall in and relax through gentle 
phase mixing. The results of our analysis in § ^ agree with this simplified description. To be more 
specific, we have found that the empirical 6c{Mo) and Xs{Mo) correlations are consistent with the 
fact that halos show, at formation, the same density profile with universal values of 6c and Xs- Until 
a new major merger takes place, the density profile of halos keeps essentially the same form, though 
the values of 5c and Xs shift as the dimensional characteristic density and scale radius, pc and r^, 



remain fixed while the virial radius R expands accordingly to the decrease of the cosmic critical 
density. As shown for the SCDM case by Avila- Reese, Firmani, & Hernandez (1997), the latter 
evolution seems to be a natural consequence of adiabatic-invariant secondary infall. On the other 
hand, some effect along the lines proposed by Syer & White (1997) might explain the universal 
halo density profiles resulting from major mergers. 

We thank J.F. Navarro for kindly providing the data for the empirical correlations appearing 
in Figures 1 and 3. The present work has been supported by the Direccion General de Investigacion 
Cientifica y Tecnica under contract PB96-0173. 
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Fig. 1. — Predicted 6c{Mo) correlations compared with the empirical data from NFW's N-body 
simulations (filled circles). Dotted and short-dashed curves show the predictions for two extreme 
values of Am, while the solid curves correspond to the value of this parameter that gives the best 
overall fit. Cosmogonies with fig < 1 contain a fourth long-dashed curve which shows, for Am = 0.6, 
the predictions arising from the assumption that pc is proportional to pcTit[z{{Mo)], instead of to 
the mean density of the universe at halo formation. 
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Fig. 2. — Distributions of log[(5c(Mo)] in the SCDM cosmogony for halos of different masses pre- 
dicted by the modified LC model with Am = 0.6. 
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Fig. 3. — Predicted Xs{Mq) correlations compared with the empirical data from NFW's N-body 
simulations (solid circles). 
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Table 1. Parameters of the models 



P{k) 


no 


Ao 


^8 


M,/Mq 


C 


4f 


a^sf^ 


a^sf 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SCDM 


1.0 


0.0 


0.63 


3.08x101''^ 


1.21x10"^ 


1.21x10"^ 


0.173 


0.229 


ACDM 


0.25 


0.75 


1.3 


6.31x10^3 


4.21x10^ 


3.77x10^ 


0.291 


0.285 


n = -1.5 


1.0 


0.0 


1.0 


1.47x10^^ 


8.30x10^ 


8.30x10^ 


0.204 


0.223 


n = -1.0 


1.0 


0.0 


1.0 


2.48x10^^ 


1.28x10"^ 


1.28x10"^ 


0.169 


0.181 




0.1 


0.0 


1.0 


2.82x10^3 


2.65x10"^ 


1.00x10"^ 


0.188 


0.184 


n = -0.5 


1.0 


0.0 


1.0 


3.40x10^^ 


2.16x10^ 


2.16x10^ 


0.135 


0.148 


n = 0.0 


1.0 


0.0 


1.0 


4.19x10^^ 


6.19x10^ 


6.19x10^ 


0.088 


0.096 




0.1 


0.0 


1.0 


4.56x10^3 


5.77x10^ 


1.33x10^ 


0.064 


0.065 


^implied by 5ci 

















